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Abstract 

We analyse the B — » Kir decays using the factorisation model with final state 
interaction phase shift included. We find that factorisation seems to describe 
qualitatively the latest CLEO data. For a test of the factorisation model, we 
derive a relation for the branching ratios independent of the strength of the 
strong penguin interactions. This relation gives a central value of (0.60 x 10 -5 ) 
for B(B° — ► K°ir°), somewhat smaller than the latest CLEO measurement, 
but the experimental errors are yet too big to take it as a real prediction of the 
factorisation model. We also find that a ratio obtained from the CP-averaged 
B — > Kir decay rates could be used to test the factorisation model and to 
determine the weak angle 7 with more precise data, though the latest CLEO 
data seem to favor 7 in the range of (90° — 120°). 
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Typeset using REVTpjX 



One of the possibilities offered by the B — > Ktt decays is the determination of the CP- 
violating phase 7, one of the angles in the (db) unitary triangle of the Cabibbo-Kobayashi- 
Maskawa (CKM) quark mixing matrix in the standard model Infact the large CP- 
averaged branching ratio(£>) for B — > Ktt as observed by the CLEO Collaboration indi- 
cate that the penguin interactions contribute a major part to the decay rates and provide 
an interference between the Cabibbo-suppressed tree and penguin contribution resulting in 
a CP-asymmetry between the B — > Ktt and its charge conjugate mode. The CP-averaged 
decay rates depend also on the weak phase 7 and give us a determination of this phase once 
a reliable description of the B — > Ktt decays could be established [f|J|. 

With the latest measurement by the CLEO collaboration , we have now the CP- 
averaged branching ratios for all the B — > Ktt decay modes. In particular, the B° — > K°tt° 
mode is found to have a large branching ratio of (l^Gt^l^) x 10~ 5 compared with a value 
in the range (0.5 — 0.74) x 10~ 5 in the factorisation model The predicted values for other 
modes are, however, more or less in agreement with experiment. As the effective Hamilto- 
nian for B — > Ktt decays is well established with the short-distance Wilson coefficients for 
tree and penguin operators now given at the next-to-leading logarithms(NLL) QCD radia- 
tive corrections the most important theoretical uncertainties would probably come 
from long-distance matrix elements obtained with the factorisation model and final state 
interaction (FSI) effects. Infact one of the main uncertainties in the penguin contributions 
to B — > Ktt decays come from the value of the current s quark mass which is not known to 
a good accuracy. There are also non-factorisation terms which must be included in the form 
of an effective Wilson coefficients to make the amplitudes scale-independent ]7],[13[]. Thus 
a more precise test of factorisation is to consider quantities which are independent of the 
strong penguin contributions. This is the main purpose of this paper. When all the B — ► Ktt 
decay modes are measured with good accuracy, and if the rescattering phase is known the 
dominant strong penguin contribution could be determined from the measured branching 
ratios assuming factorisation for the small tree-level and electroweak penguin terms, as will 
be discussed in the following. Though the present data are not yet sufficiently accurate for a 
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determination of the effective Wilson coefficients in B — > Kir decays at this time, a first step 
toward an understanding of B — > Kn decays is to see how well these penguin-dominated 
charmless B decays can be described by factorisation using the Wilson coefficients obtained 
from perturbative QCD. As argued in ||14|| , for these very energetic decays, because of color 



transparency, factorisation should be a good approximation for B — > Kir decays if the Wil- 
son coefficients are evaluated at a scale /z = 0(mb). We could thus proceed to the test of 
factorisation bearing in mind that there are possible scale-dependent corrections from non- 
factorisation terms to be determined with more precise data. To include FS1 effects, as in 
0, we assume that elastic FS1 effects can be absorbed into the two AI = 1/2 and AI = 3/2 
elastic icK — > icK rescattering phases 5i and 63 taken as free parameters and include only 
inelastic effects coming from the charm and charmless intermediate state contributions to 
the absorptive part of the decay amplitudes. These inelastic contributions can be included 
in the Wilson coefficients of the penguin operators which now have an absorptive part and 
are given in [jT0|,|12|,|T5]. 



We begin by first giving predictions in factorisation model for the B — > Kir decay rates 
and branching ratios in terms of the rescattering phase difference 5 and for a typical value 
of the weak phase 7. As will be seen, factorisation seems to produce sufficient B — > Kir 
decay rates. We could thus proceed to a test of the factorization model by comparing 
with experiments, quantities obtained by factorisation which are independent of the strong 



rescattering phase difference ||16|| . We find that the sum of the CP-averaged branching ratios 
B(B~ -> K-tt°) + B(B- -> K°n-) and B(B° -> K-ir + ) + B(B° -> AV) are independent of 
the FSI rescattering phase. Other quantities obtained from various combination of the decay 
rates, for example, the quantity A defined as T(B~ — > K ir~) + T(B° — > K~ir + ) — 2(T(B~ — > 
A~7T°) + T(B° — > A°7r )) is independent of the strong penguin contributions and could be 
used to predict B(B° — ► A°7r°) in terms of the other measured branching ratios. As the 
main purpose of this paper is to test the factorisation model using relations independent 
of the strong penguin interactions, we will not discuss here a recent theoretical work on 
factorisation in B — ► tttt decays which should be completed to have all the logarithms of rrib 



under control [|17J . 

In the standard model, the effective Hamiltonian for B — ► Kn decays are given by 
HQ 

H cS = 9l=[V ub v: s (ciO? + c 2 O u 2 ) + VaVZfaOt + c 2 O c 2 ) 
10 

-E(W c *)°*l+ h - c -. C 1 ) 

i=3 

in standard notation. At next-to-leading logarithms, q take the form of an effective Wilson 
coefficients c° s which contain also the penguin contribution from the c quark loop and are 
given in 
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The tree level operators 0\ and 2 as well as the electroweak penguin operators Oj — 0\q 
have both 1 = and 1 = 1 parts while the QCD strong penguin operators O3 — 0$ have 
only 1 = terms. The B — > Kit decay amplitudes can now be expressed in terms of the 
decay amplitudes into 1=1/2 and 1 = 3/2 final states as 0, 

Ak-^o = ^B 3 e iS * + + B 1 )e iS \ 

Ago*- = ^B 3 e^ - s^{A x + BOe* 

A*-** = ^W* 3 + ^(A 1 - B 1 )e i5 \ 

A ROn0 = 2 -B z e^ - ^(A 1 - B x )e^ , (2) 

where Ai is the sum of the strong penguin Af and the 1 = tree level Aj as well as the 
1 = electroweak penguin A^ contributions to the B — > i^7r 1=1/2 amplitude; similarly 
-Bi is the sum of the J = 1 tree level Bf and electroweak penguin contribution to the 
1 = 1/2 amplitude, and B 3 is the sum of the 1 = 1 tree level Bj and electroweak penguin 
BY contribution to the 1 = 3/2 amplitude. 

The factorisation approximation is obtained by neglecting in the Hamiltonian terms 
which are the product of two color-octet operators after Fierz reordering of the quark fields. 
The effective Hamiltonian for non-leptonic decays are then given by Eq. (|1|) with Cj replaced 



by dj and Oj expressed in terms of hadronic field operators. In the notation of Ref . || , we 
have 

^1 = ^ VubV* s r a 2 , 



'< — 7= V ub V* r 



— -a 2 + a x X 



Bj = 1 - V ub V* s r [a 2 + a x X] , 
/3 

A\ = -i— V tb V* r [a 4 + a 6 Y] , 

^/3 

A T = ~^ V tb V t :r[a 8 Y + a 10 ] 



Bt = #j 



B? = i^V tb V t :r 



B. 



w 



VtbVts r [(a 8 Y + a w ) - (a 7 - a 9 ) X] 



(3) 



where r = G F f K F*« (m 2 K )(m%-ml), X = (fJfK)(F BK (ml)/F ^(m 2 K ))(m%-m 2 K )/(m 2 B - 



m: 



2 ), Y = 2m 2 K /[{m s + m q ){mh — m q )] with q = u, d for n ± '° final states, respectively, and 



a.,- are given in terms of the effective Wilson coefficients c° ff (iV c is the number of effective 
colors) by 



cf + df+jN c for j = 1,3,5,7,9 



Gj = cf + cf_ l /N c for j = 2, 4, 6, 8, 10 



(4) 



In our analysis, we use N c = 3 and m b = 5.0 GeV which give aj the following numerical 
values 



ai = 0.07, a 2 = 1.05, 

a 4 = -0.043 - 0.0162, a 6 = -0.054 - 0.016i, 

a 7 = 0.00004 - 0.00009^, a 8 = 0.00033 - 0.00 032, 

a 9 = 0.00907 - 0.000092, a w = -0.0013 - 0.000032'. 



(5) 



Note that a± is sensitive to N c and is rather small for N c = 3 . As there is no evidence for a 
large positive a\ in B — > K7r decays which are penguin-dominated and are not sensitive to 



ai, we use ai evaluated with N c = 3 given in Eq.@. Indeed, the predicted branching ratios 
remain essentially unchanged with a\ = 0.20 taken from the Cabibbo-favored B decays 

In the absence of FSI rescattering phases, we recover the usual expressions for the decay 
amplitudes in the factorisation approximation. We have used c^ ff given at next-to-leading 
order in QCD radiative corrections [H,Pi|l2l and evaluated at a scale \x = m b . We note that 
the coefficients c| ff , c| ff , c| ff and Cg ff are enhanced by the internal charm quark loop due 
to the large time-like virtual gluon momentum q 2 = ml/2 as pointed out in H^, [10| , [l5| (the 
other electroweak penguin coefficients like Cy and Cg ff are not affected by this charm quark 
loop contribution in any significant amount). This enhancement of the strong penguin term 
increases the decay rates and bring the theoretical B — > Kir decay rates closer to the latest 
CLEO measurements. In the above expressions, the tree level amplitudes are suppressed 
relative to the penguin terms by the CKM factor V ub V* s /V t bV^ s which can be approximated 
by — (|Kifc|/|Kife|) x (| V^|/|Kid|) exp(— i 7) after neglecting terms of the order 0(A 5 ) in the 
(bs) unitarity triangle. The B — * Kir decay rates then depend on the FSI rescattering 
phase difference 5 = 5% — 5i and the weak phase 7 . In the following, we shall use the 
set of parameters of 0J2T| which give f w = 133 MeV, f K = 158 MeV, F B7r (0) = 0.33, 
F^ K (0) = 0.38. We use m s = 120MeV, \V cb \ = 0.0395, \V cd \ = 0.224 and \V ub \/\V cb \ = 0.08 
HI. At the moment, m s is not known to a good accuracy, but a value around 100 — 120 MeV 



inferred from m^* — m p , m D + — m D + and m#o — mgo mass differences |21[] seems not 
unreasonable. To show the factorisation predictions and the dependence of the branching 
ratios on the rescattering phase difference 5, we give, as an example, the CP-averaged 
B — > Kn decay rates in Fig.l evaluated with a CKM value given by p = 0.12, r\ = 0.34 
corresponding to 7 = 70°. 

As can be seen from Fig.l, all the B — > Kit decay modes for B~ and B° , except the 
B° — > K 7t° mode, have branching ratios more or less in agreement with the latest CLEO 
data OH which give, for the CP-averaged branching ratios 
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B(B + - 


-> K + n [> ) 


/^ 04-1 4\ .-, ^ f\ 

= {n.6±l-° 7 ±{i) x 10 6 , 




B(B + - 


-> K°TX + ) 


= (18.2+to± 1-6) x 1(T 6 , 




B(B° - 


> K + TT-) 


= (17.2i|| ±1.2) x 1(T 6 , 




B(B° - 


-> K°n°) 


= (14-61^) x 1CT 6 . 


(6) 



The computed decay rates shown above could be larger with the form factors given in 



22| and could bring the B — > Kn decay rates closer to the latest CLEO data. 

We now turn to the test of factorisation in B — > Kit decays. The decay rates into a Kn 



final state is given by 



T(B -> Kn) = C\A Kn \ 



(7) 



where the subscript Kir refers to any of the decay modes for B~ and B° and C is the usual 
phase space factor. By summing over the decay modes for B~ and for B° respectively, we 
have in terms of Ai, B\ and B 3 , 



b- 
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~2 
_3 


B 3 


2 + 


|^i + fill 2 
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B 3 


2 + 
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7T 



r(s- 



ir°7r-) and T 



r( J B c 



7T 



r^ 



where T B - = T(B 
K°7r°). 

The quantities T B - and T B o are independent of the rescattering phase difference S. They 
are given in the factorisation model as a function of the weak phase 7. Two other related 
quantities of interest obtained from the above Eq.(§|) are 

-2, 



r b B B - +B B o = 2C 



\B 3 \ 2 + |Ax| 2 + \B±\ 2 



T B 



r b B B - ~ B B o = 4 C Re {A\B{]tbo 



(9) 



which, together with one measured B — ► i^7r branching ratio, would enable us to deter- 
mine the strength of the strong penguin contribution as well as its absorptive part and 7 
assuming factorisation for the small tree-level and electroweak penguin contributions, if the 



rescattering phase difference 8 could be inferred from the ^-dependent branching ratio and 
from other sources. In the above expression, t b o is the B° lifetime and r& = t b o / t b - , 

Also, if all the four B —>■ Kit decay rates (CP-averaged) are measured with good accuracy, 
in particular with a precise measurement of the B° — > K°7r° branching ratio, the following 
quantities 

fl ' = fe-* = av^o (10) 

could also be used to test factorisation. 

As the CP-averaged B —>■ Kir decay rates depend on 7, the computed partial rates T B - 
and T B o would now lie between the upper and lower limit corresponding to cos(7) = 1 
and cos(7) = —1, respectively. As shown in Fig. 2, where the corresponding CP-averaged 
branching ratios (B B o and B B -) for T B ~ and T^o are plotted against 7, the factorisation 



model values with the BWS form factors |2(| seem somewhat smaller than the CLEO central 
values by about 10 — 20%. Also, B B - > B B o while the data gives B B - < B B o by a small 
amount which could be due to a large measured B° — > K n decay rates. We note that 
smaller values for the form factors could easily accommodate the latest CLEO measured 
values, if a smaller value for m s , e.g, in the range (80 — 100) MeV is used. What one 
learns from this analysis is that B — > decays are penguin-dominated and the strength 
of the penguin interactions as obtained by perturbative QCD, produce sufficient B — > Kir 
decay rates and that factorisation seems to work with an accuracy better than a factor of 
2, considering large uncertainties from the form factors and possible non-factorisation terms 
inherent in the factorisation model. With more precise measuremnts expected in the near 
future, it might be possible to have a detailed test of factorisation and a determination 
of 5 and 7 by comparing with experiments various relative branching ratios, to reduce 
uncertainties from form factors and CKM parameters. Other test of factorisation could also 
be done by looking for quantities which are independent of the strong penguin interactions. 
Infact, since the four B — > Kit decay rates depend on only three amplitudes A\, B\ and 
S3, it is possible to derive a relation between the decay rates independent of A\. From the 
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following quantities, 



7T 



X 



iBsr+ci^ir+iSii 



r(s _ 



2 



C 2 



x 



^|5 3 | 2 + (l^l 2 + l^l 2 ) + -l He (2J*iJ lC «) 



(11) 



where Ci 



|C and C 2 = |C, we obtain 



- 2 [r(£- -> ^-7T°) + T(B° -> i?°7T )] } T B o 

(Cr B o) 



-^3 

3 I 3 



2 - — Re(£ 3 *S ie i5 ) 



(12) 



From Eq. (|T2]) , we see that A is independent of A\ and hence is independent of the strong 
penguin term. It is given by the tree-level and electroweak contributions which are much 
smaller than the strong penguin term. As can be seen from Fig. 2, where its values for 5 = 
are plotted against 7. A is of the order O(10~ 6 ) compared with 13 b- and B B o which are 
dominated by the strong penguin contribution and are in the range (1.6 — 3.0) x 10~ 5 . As 
the variation with 5 is negligible, A remains at the O(10~ 6 ) level for other values of 5 ^ 0. 
Thus, to this level of accuracy, we can put A ~ 0. Eq.(^) becomes 



+ 13 



2 [B R07T o + r b B K ~ n o] 



(13) 



This relation can be used as a test of factorization with more precise measurements of the 
CP-averaged branching ratios. Conversely, it can also be used to predict B(B° — > K n ) in 
terms of the other measured branching ratios. From the latest CLEO data, with A ~ 0, 
Eq.(|l3"D then gives B(B° — ► /? 7r°) = 0.60 x 10 -5 . As can be seen, the large experimental 
errors prevent us from drawing any firm conclusion on the validity of factorisarion, though, 
the above predicted central value for B(B° — > K°7T°) is somewhat smaller than the CLEO 
data. 
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For another test of factorisation and a determination of 7 we have derived a relation in 

the form of a ratio which is independent of the form factors and the CKM parameters. It is 

given by the ratio R of the two CP-averaged quantities as 

[B(B~ -> K-tt°) + B(B~ -> K%-)] 
R = B(B- -> KH-) + B(B° -> K-n+)/r b ' ^ 

Numerically, we find that terms proportional to cos(5) and sin(5) in R is of the order 10 -7 

and thus can be safely ignored. Thus R is a function of 7 alone and can be used to determine 

7 as it does not suffer from the uncertainties in the form factors and in the CKM parameters. 

In Fig. 3 we give a plot of R as a function of 7. As can be seen, it is not possible to deduce a 

value for 7 with the CLEO data which gives R = (0.80 ± 0.25) as the theoretical prediction 

for R lies within the experimental errors. If we could reduce the experimental uncertainties 

to a level of less than 10%, we might be able to give a value for 7. Thus it is important 

to measure B — > Kit decays branching ratios to a high precision. It is interesting to note 

that the central value of 0.80 for R corresponds to 7 = 110°, close to the value (HSl^) 

found by the CLEO Collaboration in an analysis of all known charmless two-body B decays 

with the factorisation model 0] . It seems that the CLEO data favors a large 7 in the range 

(90° — 120°). A large 7 as shown in Fig. 2, would increase the factorisation values for Bb- 

and Bbo which are given numerically by 

B B - = (2.757- 0.409 cos(7)) x 10~ 5 
B B o = (2.270 - 0.624 cos(7)) x 10~ 5 

For the ratio R, we have 

R _ (2.651 - 0.393 cos(7)) 
~ (3.253 - 0.652 cos(7)) 

Also shown in Fig. 3 are the ratio R\ and R2 defined in Eq. ([TOP . As R\ shows strong 

dependence on 7, a better way to determine 7 would be to use R\ rather than R when a 

precise value for B(B° — > K n ) will be available. 

Given 7 = 110°, all the B — > Kit branching ratios can be predicted in terms of the 

rescattering phase difference 5 as shown in Fig. 4. Comparing with Fig.l, we see that, 
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(15) 
(16) 



except for B(B° — > A'°7r°) which remains at the 6.5 x 10~ 6 level, the other branching ratios 
become larger with 7 = 110° and closer to the CLEO data which indicate B~ — > K°7T~ 
and B° — > K~ir + are the two largest modes with near-equal branching ratios in qualitative 
agreement with factorisation. Fig.l shows that these two largest branching ratios are quite 
apart, except for 5 < 50° while Fig. 4 suggests S should be large, in the range of (40° — 70°). 
With a smaller 7 < 110° and some adjustment of form factors, the current s quark mass and 
CKM parameters, it might be possible to accommodate these two largest branching ratios 
with a smaller 5. We note that the dependence of the four branching ratios shown in Fig.l 
and Fig.4 are essentially the same and is given by (4/3v^3) Re(AiBg exp(i5)), apart from the 
sign, as the interference term Re(BiBg exp(i5)) is much smaller than Re(AiB:!; exp(i5)). 

We note that we have also considered a possible contribution from inelastic rescatter- 
ing effects as an additional small absorptive contribution Aj to A\ from D D* s and other 
intermediate states to the S'-matrix unitarity relation. We find that the variation of 
and T B o as a function of Aj is negligible. For this reason, we have set Ai = 0. Also, since 
the theoretical values for the decay rates shown above show qualitative agreement with the 
measured values, the strong penguin terms with enhancement by the internal c-quark loop 
seem to produce sufficient decay rates, a large dispersive inelastic contribution would not be 
needed in B — > decays. 

The CP-asymmetries, plotted against 5 as shown in Fig.5, are given by 

As CP = ~ (17) 

where T is the decay rate. The predicted CP-asymmetry Asqp for the B — > Kir decay 
modes are in the range ±(0.04) to ±(0.3) for the preferred values of 5 in the range (40 — 70)° 
mentioned above. The latest CLEO measurements [^] however, do not show any large 
CP-asymmetry in B — > Kit decays, but the errors are still too large to draw any conclusion 
at the moment. 

In conclusion, factorisation with enhancement of the strong penguin contribution seems 
to describe qualitatively the B — > Kir decays, although the predicted B(B° — > K°ti ) is 
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below the measured value. Further measurements will enable us to have a more precise test 
of factorisation and a determination of the weak angle 7 from the FSI phase-independent 
relations we shown above. 
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FIG. 1. B(B -» 
branching ratios B 



(b) 




Kit) vs. 5 for 7 = 70°. The curves (a), (b), (c), (d) are for the CP-avera^ 
" -» i^-vr , tfV" and B° -> i^-vr + , tfV, respectively. 
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FIG. 2. B B _ (a), B §0 (b), A (c) vs. 7 
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FIG. 3. The curves (a), (b), (c) are for R, R2, Rl respsectively 
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